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^^ ABSTRACT 

o 

vN Context. This paper is part of a larger project in which we want to focus on the still poorly understood AGB third dredge-up processes 

^ and associated s-process nucleosynthesis. 

^ Cy Aims. We confront accurate spectral abundance analyses of post- AGB stars in both the LMC and SMC, to state-of-the-art AGB model 

^I predictions. With this comparison we aim at improving our understanding of the 3rd dredge-up phenomena and their dependencies 

^^ on initial mass and metallicity 

^^ Methods. Because of the well constrained distance with respect to Galactic post- AGB stars, we choose an extra-galactic post- AGB 

^ star for this contribution, namely the only known 21 yum object of the Small Magellanic Cloud (SMC): J004441. 04-7321 36.4 . We 
used optical UVES spectra to perform an accurate spectral abundance analysis. With photometric data of multiple catalogues we 

r~i construct a spectral energy distribution and perform a variability analysis. The results are then compared to predictions of tailored 

P^ theoretical chemical AGB evolutionary models for which we used two evolution codes. 

^O Results. Spectral abundance results reveal J00444 1.04-732 136.4 to be one of the most s-process enriched objects found up to date, 

• while the photospheric C/0 ratio of 1 .9 + 0.7, shows the star is only modestly C-rich. J004441 .04-732136.4 also displays a low [Fe/H] 

fC = -1.34 + 0.32, which is significantly lower than the mean metallicity of the SMC. From the SED, a luminosity of 7600 + 200Lo is 

Mh found, together with E(B-V) = 0.64 + 0.02. According to evolutionary post- AGB tracks, the initial mass should be ^ 1.3Mo. The 

^ photometric variability shows a clear period of 97.6 + 0.3 days. The detected C/0 as well as the high s-process overabundances (e.g. 

^ [Y/Fe] = 2.15, [La/Fe] = 2.84) are hard to reconcile with the predictions. The chemical models also predict a high Pb abundance, 

^*T^ which is not compatible with the detected spectrum, and a very high ^^C/^^C, which is not yet constrained by observations. The 

^ predictions are only marginally dependent on the evolution codes used. 

i«_i Conclusions. By virtue of their spectral types, favourable bolometric corrections as well as their constrained distances, post- AGB 

stars in external galaxies off'er unprecedented tests to AGB nucleosynthesis and dredge-up predictions. We focus here on one object 

^H J004441. 04-732136. 4, which is the only known 21 jim source of the SMC. We show that our theoretical predictions match the s- 

^ process distribution, but fail in reproducing the detected high overabundances and predict a high Pb abundance which is not detected. 

CO Additionally, there remain serious problems in explaining the observed pulsational properties of this source. 

_j. Key words. Stars: AGB and post- AGB - Stars: abundances - Stars: evolution - Stars: oscillations - Galaxies: Magellanic Clouds - 

_i. Nuclear reactions, nucleosynthesis, abundances 
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p^ 1, Introduction internal chemical ev olution of these stars ( Habing & Qlofsson 

[20031 |Herwig]2QQ5] ). 



The final evolution of low- and intermediate-mass stars is a fast 



>< 



r ^u \ ^^r^-^o u/Ar^o\ ^u Here we focus on the the poorly understood AGB 3rd 

transition from the Asymptotic Giant Branch (AGB) over the , , , , . , . , , . r ., • . t 

^ A r^o ^ V ^ ^1 fu r>T ^ TVT u 1 r>u /r>TVT\ u drcdgc-uD phenomenon, during which products of the internal 

post- AGB transit towards the Planetary Nebula Phase (FN), be- ,^ ^J . i , . . ., r r ., . rr., • • 



r ^u 4^ n ^ T ^ ^T7uv TA ^ /^T7TA\ nuclcosynthcsis are brought to the surface of the star. This is 

^ fore the stellar remnant cools down as a White Dwarf (WD). . i i?^ , . , ^ , , , ^ 

?-H Ai^u u ^u- u u IT 1 T ^ ^ ^u mainly '^C as the primary product or the triple alpha reaction, 

rri Although this scheme may be generally acknowledged, there , ^ / ^, i : r . . .1-4,, ^ 

^^. j^j- ly n ^ ' ' , r -.'rr ^' ^ ^ but also thc products of ncutron-capturc synthcsis. Thcrc arc two 

IS no understanding from first principles of dmerent important . ^ . .^^ .x , ?7^t / ^?^i,^ 

, . 1 ^, ^ \u ^ 4^' u T^u mam neutron sources m AGB stars: 1) the ^^Ne(a,n)^-'Mg reac- 

physical processes that govern these evolutionary phases. The . , . , . . , ^^ ^ o\r^ i ^6?^ i 

. ,^ . i^j^^uii r ^ 4^ A' f tion which IS activated at temperatures or T > 300 x 10"K, and 

mam shortcomings are related to the lack or understanding or ^. . n^. . i^^ . , . , . .• 7 i . i i 

^, 1 u ' A 1 1 ^' T ^u 2) the ^^C((2,n)^"0 reaction, which is activated at much lower 

the mass-loss mechanisms and mass-loss evolution along the ^ ^ ' ^, ^ „ „ .^^6T^ r^^ • i i i - ^ 

Ar^o ^^uu ^u- ^^u- temperatures or T > 90 x 10"K. Observational and theoretical 

AGB ascent, the subsequent shaping processes of the circum- . f , , , , n^, .i6„ . . , 

.11 ,11 a4^u\\^^a\\ a ^ a- ^^u evidence has shown that the ^C(Qf,n)'"0 reaction is the mam 
stellar shells, and the lack or fundamental understanding or the ^ . , a^t^\ r . ^.^i, .c^. • n 

neutron source m low-mass AGB stars of ^ 1 - 3Mo (Straniero| 

TT"; ^TT- ^ ^ T^ T^ c ^. 1 let al.|1995l|Gallino et al.|1998l|Abia et al.|2002|). 

Send offprint requests to: K. De Smedt, ken- ' < " ■ f' < » 

neth.desmedt@sterkuleuven.be Synthesis by the s-process in AGB stars is an important con- 

* based on observations collected with the Very Large Telescope at tributor to the cosmic abundances past the iron peak and these 

the ESO Paranal Observatory (Chih) of programme number 084.D- stars are also thought to be very importa nt contributors to the 



0932. total carbon and nitrogen enrichment (e.g. Romano et al.|201(^ 
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Kobayashi et al.|20TT] ). Post-AGB photospheres bear witness to 
the total chemical changes accumulated during the stellar life- 
time. 

In recent years, theoretical models of these internal nucle- 
osynthesis and photospheric enrichment processes gained enor- 
mously in sophistication: extensive nuclear networks with up- 
dated cross-sections were included (e.g. [Cristallo et al. 2Q11[ 
[20091 [Karakas.2010, Church et al.|[2nDgl ); non-convective mix- 
ing such as differential rotation and thermohaline mixing have 
been imp lement ed (e.g.|Siess et al.|2003t |Siess|2007t |Stancliffe| 
et al.|200 7 ; Ange lou et al.|201 1| ); the effect of deep mixi ng or ex- 



tra m ixing proces ses have bee n critically evaluated (e.g. Karakas 
|et al.,2010^ Busso et al.|2010| ) and overshoot regimes have been 
explore d in more deta il in order to explain the extent of the ^^C 
pocket ( |Herwig|2005 ). All these processes involve great uncer- 
tainties and observational data are required to calibrate them. 

Accurate determi nation of photosph eric abundances in AGB 
stars is difficult (e.g. |Abia et ari|2008| and references therein). 
First of all, the photosphere is dominated by molecular opacity, 
making the mere detection of trace elements over a wide range 
of atomic masses difficult. Moreover, AGB stars often have dy- 
namic atmospheres caused by pulsations which develop into dust 
driven winds. Finally, chemical results of AGB stars are difficult 
to interpret, because stars with different initial mass and metal- 
licity occupy the same region in the HR-diagram. 

Post-AGB stars do not have these drawbacks: First, their at- 
mospheres do not show the large amplitude pulsations as well as 
the large mass-loss rates that characterise evolved AGB atmo- 
spheres. Second, their photospheres are hotter, so atomic tran- 
sitions prevail. This allows to quantify the abundances in post- 
AGB stars for a very wide range of elements, from CNO up to 
the most heavy s-process elements, well beyond the Ba peak 



(e.g. I Van Winckel & Reyniers||200()| Reyniers & Van Winckel 
|2Q03| r 

During the past decade it has been realised that Galactic post- 
AGB stars are chemically much more diverse than anticipated 
(e.g. Van Winckel||2003 ). Some post- AGB stars are i ndeed the 
most s-pro cess enriched objects known to date (e.g. Reyniers 
|et al.|2004| ), while others are not enriched at all. A distinct sub- 
class of Galactic post- AGB stars is formed by the so-called 21 
fim objects display ing a strong solid-state feature around 21 yum 
( [Kwoket al.|1989| ) in their IR spectra. The carrier of this feature 
still needs to be identified, alth ough severa l suggestions have 
been discussed in literature (e.g.|Posch et al._2004j. The feature 



is only detected in post-AGB Carbon stars (e.g. Hrivnak et al. 
|2008[ 2009| ) and till now not yet in spect ra of normal carbon stars 
nor in carbon-rich planetary nebula (Vol ketal.|2011| ). The chem- 
ical studies of Galactic 21 yum stars show that they display strong 
overabundances of s-pro cess elements with a wide range of neu- 
tron capture efficiencies (fVan Winckel & Reyniers |2000[ Reddy 



et al.|2002|[Reyniers & Van Winckel|2003t [Reyniers et al.|2004| ). 
The poorly known distances and hence luminosities and masses 
of the limited Galactic post-AGB sample hamper, however, the 
interpretation of the variety of abundances in the broader theo- 
retical context of stellar (chemical) evolution. 

We therefore initiated a project to exploit our newly iden- 
tified l arge sample of post-A GB stars in the Large Magellanic 
Cloud ([van Aarle et al.|201 1[ ) and in the Small Magellanic Cloud 
( Kamath et aLpOlf ), to study the s-process production and as- 
sociated 3rd-dredge-up processes. The unique spectral charac- 
teristics of post-AGB stars, together with the new large sam- 
ple, covering a wide range in luminosities and metallicities, with 
well constrained distances, means that these objects provide un- 
precedented direct tests for the theoretical structure and enrich- 



ment models of solar-mass stars. In this contribution, we fo- 
cus on the object J004441. 04-732136.4 (hereafter abbreviated 
to J004441.04) which is the only object of the Small Magellanic 
Cloud (SMC) known to date in which the 21yum feature is iden- 
tified fVo lk et aT]|2011 ). The paper is organised as follows: In 
Sect. [2] we discuss the observational data of J004441.04. Sect, 
^describes the spectral analyses performed to determine the at- 
mospheric parameters as well as the accurate photospheric com- 
position. The spectral energy distribution (SED) of J004441.04 
is analysed in Sect. [4] followed by the determination of the ini- 
tial mass in Sect. 5. We then used two state-of-the-art evolution 
codes coupled to post-processing codes to compute the predicted 
chemical signature of this object in Sect. [6] The discussion and 
conclusions are given in Sect. [7] 



2. Data 

2.1. Photometric data 

Photometric data from the Magellanic Cloud Photometric 
Survey (MCPS, [ZSitsky et al.|([2002^) the Deep Near-Infrared 
Survey (DENIST lFouque et al.| ( [20001), 2MASS (|CuMet~aL] 
2003) and the Spitzer S^MC surve7(Bolatt o et aL]|2007| ) are 
used for the construction of the SED of J004441.04 in Sect. [H 
A summary of the absolute photometric datapoints is given in 
Table [1] The MCPS, 2MASS and DENIS data were taken in a 
timespan of 2 years (1997 till 1999) while the Spitzer data were 
acquired in 2004 and 20 05. For t he variability anal ysis, we used 
the OGLE II lightcurve ( |Szymanski.2005, Udalski et Sil.fTWi) 
obtained in the I filter. 



2.2. Spectroscopic data 

We use h igh-resolution spectra obtained with the UVES spec- 
trograph (IDekker et al. 12000), which is the echelle spectrograph 
mounted on the 8m UT2 Kueyen Telescope of the VLT array 
at the Paranal Observatory of ESO in Chili. Multiple spectra of 
J004441.04 were obtained on the same day within a time span 
of approximately two hours. We used the dichroic beam-splitter 
resulting in a wavelength coverage from approximately 3280 to 
4530 A for the blue arm of UVES, and from approximately 4780 
to 5770 A and from 5800 to 6810 A for the lower and upper part 
of the mosaic CCD chip respectively. Each wavelength range 
was observed three times with an exposure time of 2871 seconds 
each. 

The UVES reduction pipeline was used for the reduction of 
the spectra. This includes the standard steps of extracting frames, 
determining wavelength calibration and applying this scale to 
flat-field divided data. Also cosmic clipping was included in the 
reduction. 

All spectra were normalised by fitting in small spectral win- 
dows, fifth order polynomials through interactively defined con- 
tinuum points. Once all three subspectra of the same wavelength 
range are normalised, the weighted mean subspectrum was de- 
termined. These mean subspectra are then again merged into one 
normalised spectrum and this is used for the spectral abundance 
determination. 

Unfortunately, the S/N is too poor for a large part of the 
blue spectrum (wavelength range 3280 to 4200 A) making these 
wavelength ranges unusable for an accurate spectral abundance 
analysis hence they are not used for the study of J004441.04. 
Except for these bluer wavelength ranges, an overall S/N around 
100 is obtained. 
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Fig. 1. Not-normalised UVES spectrum of J004441.04, centered 
around the Ba II line at 6496.897 A. The spectrum is set to rest 
velocity. 

Table 1. Photometric data used for the SED of J004441.04. A 
standard error of 0.05 mag is assumed for the SED study. IRAC 
and MIPS are instruments of the Spitzer satellite. 



band 


survey 


mag 


Fv (Jy) 


V 


MCPS 


15.96 


1.573 xlO-^ 


I 


MCPS 


14.46 


5.635 xlO-3 


I. 


DENIS 


14.44 


4.260 xlO-3 


J 


DENIS 


13.62 


6.310x10-3 


K 


DENIS 


12.94 


4.194x10-3 


J 


2MASS 


13.63 


5.635 xlO-3 


H 


2MASS 


13.23 


5.216x10-3 


K5 


2MASS 


13.02 


4.128x10-3 


[3.6] 


IRAC 


12.33 


3.976 xlO-3 


[4.5] 


IRAC 


11.90 


3.434 xlO-3 


[5.8] 


IRAC 


10.39 


1.010x10-2 


[8.0] 


IRAC 


8.344 


3.459 xlO-2 


[24] 


MIPS 


5.141 


5.998 xlO-2 



3. Spectral analyses 

It is useful to compare the spectrum of J004441.04 with Galactic 
post-AGB stars. In Eig.[2]and[3]we compare a small spectral part 
of J004441.04 with the same spectral range of two Galactic post- 
AGB stars. The lower spectrum is from the Galactic post-AGB 
star IRAS06530-0213 (hereafter abbreviated to IRAS06530). 
This is a post-AGB star with a metallicity of [Ee/H] = -0.46 
and Teff = 7250 K and it is recognised to be the star with the 
highest s-process overabundances known to date ( Reyniers et al. 
2004| ). The upper spectrum is from the non-enriched Galactic 
star H D 112374 = HR4912 (Lambert et al.|1983t|Giridhar et al.| 
|1997 |). This star has a similar spectral type and a metallicity of 
[Ee/H] =-1.1. The normalised spectra are brought to the same 
velocity scale and offset for clarity. Red and black vertical lines 
mark positions of s-nuclei and non s-nuclei respectively. It is 
remarkable that the spectra of J004441.04 are very similar to 
IRAS06530-0213 which illustrates its high s-process overabun- 
dance. We therefore choose the used line list of IRAS06530 in 
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Fig. 2. Comparison of the normalised spectra of J004441.04 
(middle), IRAS06530 (lower) and HDl 12374 (upper). The up- 
per and lower spectra have been shifted in flux for clarity. Red 
and black vertical lines mark positions of s-nuclei and non s- 
nuclei respectively. Eor more information, see text. 
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Reyniers et al. (2004) for quantified relative spectral analyses. 



Fig. 3. Comparison of the normalised spectra of J004441.04 
(middle), IRAS06530 (lower) and HDl 12374 (upper). The up- 
per and lower spectra have been shifted in flux for clarity. Red 
and black vertical lines mark positions of s-nuclei and non s- 
nuclei respectively. 



3.1. General methods 

Both the atmospheric parameter determination and the abun- 
dance determination are performed using the same routines. The 
local thermal equilibrium (LTE) Kurucz-Castelli atmosphere 
models ( [Castelli & Kurucz|2Q04| ) are used in combination with 
the LTE ab undance calculation routine MOOG (version July 
2009) by C. |Sneden| ( p773l ). Non-LTE efl'ects are not taken into 
account for the spectral analyses of J004441.04. 

The equivalent width (EW) of lines are measured via direct 
integration and abundances are computed by an iterative process 
in which the theoretical EWs of single lines are computed for 
a given abundance and matched to the observed EWs. Blended 
lines are avoided as much as possible. To check whether lines 
are part of blends with other identified lines, synthetic spectra 
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Table 2. Determined atmospheric parameters of J004441.04 
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Fig. 4. Comparison between B aimer lines of J004441.04 with 
synthetic models of'Coel ho et al.| POO S') with different log g. 
The upper panel shows Hq,, the lower panel shows H^. The black 
spectrum represents J004441.04, the red and green spectrum are 
synthetic spectra with a log g of 0.5 and 5.0 dex respectively. For 
more information, see text. 



are modelled with MOOG using VALD linelists ( |Kupka et al.| 
[1999 ) and compared to the stellar spectra. 

In order to identify lines, the radial velocity of J004441.04 
is estimated by fitting a Gaussian curve to a number of identi- 
fied atomic lines to determine their central wavelength. Using 
the equation of the Doppler shift, this results in a heliocentric 
radial velocity v = 148 ± 3 km/s which is accurate enough for 
line identifica tion purposes. The useful line identification tool of 
|Lobel| (|2006') was used for the first identifications. The heliocen- 
tric radial velocity of the SMC is in average 160 km/s ( Richter_ 
|et al.||1987| ) and the velocity of J004441.04 confirms the mem- 
bership of the SMC. 

3.2. Atmospheric parameter determination 

The atmospheric parameters are derived using Fe I and Fe II 
lines and the standard spectroscopic methods: the eff'ective tem- 
perature Teff is derived by imposing the iron abundance, derived 
from individual Fe I lines, to be independent of lower excitation 
potential; surface gravity log g is derived by imposing ionization 
equilibrium between the iron abundance of individual Fe I and 
Fe II; the microturbulent velocity ^t is derived by imposing the 
iron abundance from individual Fe I lines to be independent of 
reduced equivalent width (RW). For these analysis, 17 Fe I and 
4 Fe II lines are used. Due to the high s-process enrichment of 
J004441.04, the number of useful non-blended Fe lines for the 
atmospheric determination is limited, especially for Fe II. 

The results for the diff'erent atmospheric parameters are 
shown in Table [2] and lie in the range of typical post- AGE pa- 
rameter values. The iron abundance for the determined atmo- 
spheric parameters is a good estimate of the overall chemical 
composition [M/H] which is in this case [Fe/H] = -1.34 ± 0.32 
(dex). The abundance error of ± 0.32 is the total error includ- 
ing the line to line scatter and the uncertainties by the used at- 



mospheric model which will be further discussed in Sect. 3.3 



The low iron abundance of J004441.04 with respect to the mean 
metalHcity of the SMC of [Fe/H] ^ -0.7 ( |Luck et al.||T998l ) 



Teff (K) 


6250 ± 250 


log g (dex) 


0.5 ± 0.5 


^t (km/s) 


3.5 + 0.5 


[M/H] (dex) 


-1.34 + 0.32 



classifies J004441.04 as a low-metallicity star which are gen- 



erally acknowledged as astrophysical production sites of heavy 
s-nuclei provided the third dredge-up (TDU) takes place. 

An attempt was made to determine the surface gravity of 
J004441.04 by fitting synthetic spectra to its Balmer lines Hq, 
and H^ at 6562.8 A and 4861.3 A respectively. The comparison 
between the Balmer wings of the synthetic models with diff'er- 
ent log g and J004441.04 would provide an estimate of the sur- 
face gravity of J004441.04. However, this method is not usable 
for J004441.04, since its Balmer wings display emission lines 
as shown in Fig. |4] The upper panel shows the comparison be- 
tween Ha fines, the lower panel shows H^. The black spectrum 
represents J004441.04, the red and green spectrum are synthetic 
spectra of |Coelho et al.| ( [2Q05| ) with a log g of 0.5 and 5.0 dex 
respectively. Both panels display emission in the wings of both 
Balmer lines of J004441.04 which are most probably caused by 
ongoing mass loss. Also the strong enrichment of J004441.04 
makes it impossible to fit the Balmer wings, since J004441.04 
displays multiple strong atomic lines which are not included in 
the synthetic spectra. This is especially clear in the lower panel 
ofFig.g 

The atmospheric parameters can be checked using results on 
other species only when the number of lines for that species is 
significant. Unfortunately, most species only have two, three or 
four useful lines (see already Table |3]) due to the high number 
of blends. Fe I is the only neutral ion with single lines covering 
a large range in excitation potential. Fe is also the only element 
for which a useful number of lines of diff'erent ions is found. 
Therefore, the microturbulent velocity ^t is the only parameter 
which can be checked using another species than Fe: although 
the used number of lines of La II, Ce II and Nd II is small, they 
do provide a mean to check the derived ^f This ^t check with 
La II and Cell yields ^t = 3.5 km/s and Ndll gives ^t = 3.0 
km/s , which confirms the found microturbulent velocity. 

3.3. Abundance determination 

With our preferred model atmosphere as basis, we started with a 
full but strictly relative abundance analysis and limited ourselves 
first to the lines used in the spectral analyses of IRAS06530- 
0213. 

We used mainly isolated non-blended, non-saturated lines 
and the individual atomic lines are all double-checked via a spec- 
trum synthesis to investigate the possible presence of unresolved 
blends. N lines with EWs larger than 3 A are not found and are 
not used due to possible confusion with noise in the spectrum. A 
redder spectrum will be needed for the N abundance determina- 
tion. 

Unfortunately, at these high overabundances, all Sr and Ba 
lines are heavily saturated making accurate abundance determi- 
nation of these two most famous s-process species impossible 
(Figs. [T] and [3]). The other s-process abundances come from iso- 
lated single lines except for those species where all detectable 
lines turned out to be blended. These specific blends were fitted 
by creating synthetic spectra in MOOG using VALD line lists of 
specific wavelength ranges. The resulting fits for Eu II and Gd II 
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are shown in Fig.[6|for the Eu II line at 6437.640 A in the upper 
panel and the Gd II line at 5733.852 A in the lower panel. The 
black spectrum is the spectrum of J004441.04, the colored spec- 
tra are synthetic spectra with different abundances of the studied 
element. The Eu II line at 6437.640 A is part of a blend with a 
very weak Y II line at 6437.169 A, the Gd II line is part of a 
blend with a relatively strong Ce II line at 5733.692 A. Fig. J5| 
shows the spectrum synthesis of the W II line at 5104.432 A. 
The black spectrum is J004441.04, the colored spectra are syn- 
thetic spectra with different W abundances. This W II line forms 
a blend with a Sm II line at 5 104.479 A which is indicated in the 
figure. The positions of lines unknown to VALD are indicated 
with '?' proving that there are still a large number of lines that 
need to be identified in these highly enriched object. 

The next step in the process was to constrain abundances for 
elements which have no isolated single lines. We used VALD 
line lists to predict the equivalenth widths for all lines of a given 
element by assuming a strong overabundance. For each element 
we then choose the strongest line in our observed spectrum and 
fitted this line with synthetic spectra of MOOG, analog to the Eu 
II and Gd II abundances so with appropriate input abundances of 
all elements which we could quantify previously. Unfortunately, 
also here the strong overabundances of s-process elements cre- 
ate spectral lines which are unknown to VALD linelists making 
abundance determinations difficult because of possible blends of 
unknown origin. Using the synthetic spectrum method, we esti- 
mated abundances for Mg I and Zn I together with the s-process 
elements Dy II, Er II, Yb II, Lu II and Hf II. Afterwards, we de- 
termined the EW of each used line in order to perform an error 
analysis on the derived abundances. For the abundance deter- 
mination of O, we did not use the forbidden O I lines at 6300 
and 6363 A as they are part of an identified and unidentified 
blend respectively. After the abundance determination, we per- 
formed a spectral synthesis of the 6300 A line using the derived 
O abundance which rendered a good fit. As forbidden lines are 
not sensitive to non-LTE effects (see e.g. |Kiselman|p002| )), this 
is a confirmation of the obtained high oxygen abudance. 

The final abundance results of J004441.04 are shown in 
Table |3] where N represents the number of lines used for the 
abundance determination of the species, the full line list used for 
the abundance determination is shown i n Tab le \5\ Solar abun- 
dances are taken from [Asplund et al.| p009| ), uncertanties in 
log 6 and log 6© due to line-to-line scatter are given by cr^ 
and (Te^ respectively. The sensitivity of [X/Fe] to uncertainties 
in the line-to-line scatter, effective temperature, surface grav- 
ity and microturbulent velocity are given by cr[x/Fe] ^ 0-7^// ^ (^logg 
and cr^^ respectively. The combined uncertainty in [X/Fe] due to 
the above described sensitivity is represented by crtot. Elemental 
abundances derived via synthetic spectrum fitting are indicated 
in an italic font in Table [3] We remark again that non-LTE ef- 
fects are not taken into account in the abundance determination 
of the different species. Although the abundances of some use- 
ful elements for nucleosynthesis studies cannot be determined. 
Table [3] still contains quantified abundances of a wide range of 
s-elements. 



3.4. Abundance results 

The [X/Fe] results of Table |3] are plotted in Fig.|7j the errorbars 
represent the total uncertainty (Ttot- For clarity some elements are 
labelled. The figure illustrates well the post-carbon star signature 
of J004441.04, with a C enrichement and very strong s-process 
enrichment. The C/0 ratio is 1.9 ± 0.7. The uncertainties of the 
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Fig. 5. Spectrum synthesis of the W II line at 5104.432 A. The 
black spectrum is J004441.04, the colored spectra are synthetic 
spectra with different W abundances. The Sm II line at 5104.080 
A, which forms a blend with W II, is indicated together with 
some atomic lines which are unknown in the line list. For more 
information, see text. 



model atmosphere has the biggest impact on the accuracy of this 
ratio. J004441.04 is a post-AGB Carbon star with a C/0 larger 
than one, but only mildly so. 

Concerning the available a-elements Mg, Si, S, Ca and Ti, 
the simple mean of the [X/Fe] is [a/Fe] = +0.4. Such an 
enhancement is normal for Galactic objects in this metallicity 
range as a consequence of galactic chemical evolution and would 
not point to an intrinsic enhancement. However, it is not clear 
whether this relation also holds for the SMC. 

The extreme s-process abundances point to a very effective 
dredge-up. With [Y/Fe] = +2.2 and [La,Ce,Pr,Nd/Fe] all inbe- 
tween 2.5 and 2.8 the s-process overabundance is extreme. Also 
the s-process overabundances of elements well beyond the Ba 
peak are high. The abundance results show that J004441.04 is 
the strongest s-process enriched extra-galactic object studied to 
date. 

Fig. |7] shows a stronger overabundance of the heavy s- 
process component (magic neutron number 82) around La and 
Ce with respect to the light s-process component (magic neutron 
number 50) around Y and Zr which indicates an effective ^^C 
pocket for the production of heavy elements. 

To get a good comparative view of the s-process nucleosyn- 
thesis of J004441.04 and IRAS06530, we scale the [Zr/Fe] result 
of J004441.04 to the [Zr/Fe] result of IRAS06530 making both 
[Zr/Fe] values overlap. This difference between the [Zr/Fe] re- 
sults is then added to the ratios of J004441.04. 

The results of this scaling are shown in Fig. [8] in which the 
red diamonds are the [X/Fe] results of J004441.04, open cir- 
cles represent the results of IRAS06530. The red errorbars of 
J004441.04 ratios represent the total uncertainty in [X/Fe] while 
for IRAS06530, the bl ack errorbars represent only the line-to- 
line scatter taken from Reyniers et al.| ( |2004| ). Some elements 
are labelled for clarity. Despite a metallicity difference. Fig. [8] 
clearly shows that both objects have very similar s-process abun- 
dance patterns. With a C/0 ratio of 2.8 and 1.9 respectively, also 
the C/0 ratio of IRAS06530 and J004441.04 are very similar. 
The same abundance trends can be found for the lighter ele- 
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Table 3. Abundance results of J004441.04 together with the calcuted errors due to line to line scatter and uncertainties of the chosen 
atmosphere model. N represents the number of lines used for the abundance determination of the species. Uncertainties in log e 
and log 6© due to line-to-line scatter are given by cr^ and cr^ respectively. cr[x/Fe]^ ^T,ff, ^logg ^i^^ <^^r represent the sensitivity of 
[X/Fe] to uncertainties in the line-to-line scatter, effective temperature, surface gravity and microturbulent velocity respectively. (Jtot 
represents the t otal uncertainty of [X/Fe] due to the line-to-line scatter and model uncertainties. Solar abundances are retrieved from 
Asplund et al. ( 2009) . Abundances of elements indicated with italicised text are determined via synthetic spectra fitting. For more 
information, see text. 



Species N log 6 + o-^ log eq ± a, [X/Fe] + atot o-[x/Fe] 



0-T,ff CTlogg Cr^f 

+250 K -250 K +0.5 dex +0.5 km/s -0.5 km/s 



CI 


6 


8.76 + 0.10 


8.43 + 0.05 


1.67 + 0.36 


0.14 


-0.23 


0.24 


0.11 


0.02 


-0.02 


01 


3 


8.49 + 0.09 


8.69 + 0.05 


1.14 + 0.50 


0.13 


-0.35 


0.32 


0.12 


0.02 


-0.05 


Nal 


2 


5.85 + 0.08 


6.24 + 0.04 


0.95 + 0.14 


0.12 


-0.08 


0.06 


-0.02 


0.01 


-0.03 


Mgl 


1 


6.75 ± 0.20 


7.60 ± 0.04 


0.49 ± 0.26 


0.22 


-0.11 


0.05 


-0.05 


0.02 


-0.04 


Si I 


2 


6.63 + 0.01 


7.51+0.03 


0.46 + 0.16 


0.09 


-0.11 


0.07 


-0.05 


0.02 


-0.04 


SI 


2 


6.07 + 0.13 


7.12 + 0.03 


0.29 + 0.26 


0.16 


-0.17 


0.15 


0.03 


0.02 


-0.04 


Cal 


4 


5.33 + 0.03 


6.34 + 0.04 


0.33 + 0.06 


0.09 


-0.03 


0.03 


0.01 


0.00 


0.00 


Sell 


3 


2.18 + 0.08 


3.15 + 0.04 


0.35 + 0.14 


0.15 


0.08 


-0.08 


0.01 


-0.01 


0.01 


Till 


3 


4.02 + 0.21 


4.95 + 0.05 


0.39 + 0.16 


0.25 


0.06 


-0.05 


0.01 


-0.02 


0.01 


CrII 


3 


4.26 + 0.04 


5.64 + 0.04 


-0.06 + 0.08 


0.13 


-0.02 


0.00 


-0.02 


0.01 


-0.02 


Mnl 


2 


4.77 + 0.11 


5.43 + 0.04 


0.68 + 0.14 


0.14 


-0.07 


0.01 


-0.05 


0.01 


-0.04 


Fel 


17 


6.16 + 0.07 


7.50 + 0.04 


0.00 + 0.03 


0.11 


0.00 


0.00 


0.00 


0.00 


0.00 


Fell 


4 


6.18 + 0.11 


7.50 + 0.04 


0.00 + 0.08 


0.17 


0.00 


0.00 


0.00 


0.00 


0.00 


Nil 


4 


5.09 + 0.17 


6.22 + 0.04 


0.21+0.13 


0.19 


-0.05 


0.01 


-0.04 


0.03 


-0.04 


Znl 


1 


3.85 ± 0.20 


4.56 ± 0.05 


0.63 ± 0.23 


0.22 


-0.06 


0.00 


-0.04 


-0.01 


0.00 


YII 


2 


3.04 + 0.03 


2.21 + 0.05 


2.15 + 0.42 


0.13 


0.26 


-0.08 


0.20 


-0.13 


0.19 


Zrll 


2 


3.23 + 0.11 


2.58 + 0.04 


1.97 + 0.22 


0.17 


0.13 


-0.08 


0.07 


-0.06 


0.07 


La II 


9 


2.62 + 0.13 


1.10 + 0.04 


2.84 + 0.32 


0.18 


0.23 


-0.14 


0.09 


-0.08 


0.10 


Cell 


8 


2.79 + 0.13 


1.58 + 0.04 


2.53 + 0.18 


0.18 


0.12 


-0.12 


-0.01 


-0.02 


0.01 


PrII 


2 


2.10 + 0.01 


0.72 + 0.04 


2.70 + 0.32 


0.12 


0.23 


-0.18 


0.02 


-0.06 


0.06 


Ndll 


7 


2.84 + 0.10 


1.42 + 0.04 


2.74 + 0.31 


0.16 


0.24 


-0.17 


0.03 


-0.05 


0.07 


Smll 


2 


1.85 + 0.04 


0.96 + 0.04 


2.21+0.30 


0.13 


0.22 


-0.16 


0.05 


-0.05 


0.07 


EuII 




1.13 ±0.10 


0.52 ± 0.04 


1.93 ±0.24 


0.16 


0.10 


-0.13 


-0.06 


-0.04 


0.04 


Gdll 




1.75 ±0.10 


1.07 ±0.04 


2.00 ± 0.20 


0.16 


0.05 


-0.10 


-0.04 


0.00 


-0.02 


Dyll 




2.05 ± 0.20 


1.10 ±0.04 


2.27 ± 0.22 


0.12 


0.11 


-0.16 


-0.05 


-0.01 


-0.01 


ErII 




1.90 ±0.20 


0.92 ± 0.05 


2.30 ± 0.25 


0.24 


0.04 


-0.07 


-0.04 


0.01 


-0.02 


Ybll 




2.15 ± 0.20 


0.84 ±0.11 


2.63 ± 0.63 


0.26 


-0.02 


-0.01 


-0.04 


0.57 


0.00 


Lull 




1.30 ±0.20 


0.10 ± 0.09 


2.52 ± 0.15 


0.15 


0.02 


-0.06 


-0.03 


-0.04 


0.06 


Hfll 




2.15 ± 0.20 


0.85 ± 0.04 


2.62 ± 0.13 


0.12 


0.03 


-0.08 


-0.03 


0.00 


-0.01 


WII 




2.25 ± 0.20 


0.85 ± 0.12 


2.72 ±0.27 


0.26 


0.02 


-0.07 


-0.04 


-0.01 


-0.01 



ments where only Mn displays a strong diff'erence. Following 
the index definitions of Reyniers et al. ( 2004 ), we find an intrin- 
sic index [hs/ls] = +0.5 and [s/Fe] = 2.4. In these indices, hs 
stands for high-mass s-process elements around the Ba peak and 
Is for elements around the Sr peak. We quantified the Ba abun- 
dance using the observed Nd abundance. This places J004441.04 
in the upper right corner of the upper panel in Fig. 7 of Reyniers 
et al. ( 2004| ) confirming the correlation between the total enrich- 
ment of s-process elements and the [hs/ls] index. The results of 
J004441.04 corroborate the finding that there is no clear correla- 
tion between metallicity and the [hs/ls] index in the metallicity 
range sampled by post-AGB stars ( [Reyniers et al.|2004) . 



4. Spectral energy distribution 

The spectral energy distribution (SED) of J004441.04 and its as- 
sumed distance provide the opportunity of determining its lumi- 
nosity. The photometric data of Table [T] are used for the SED 
construction of J004441.04. For post-AGB objects in the SMC, 
we need to take into account three possible sources of reddening. 
The first one is reddening by interstellar dust in the Milky 
Way Galaxy (MW) towards the SMC. |Schlegel et al.| ( |1998| ) de- 



rived the extinction towards the SMC by galactic dust to be small 
with an average of E(B-V) = 0.037 mag. 

The second possible reddening source is reddening in the 
SMC itself. Fig. [9] shows the MW (g reen) and SMC (blue ) ex- 
tinction curves whic h are taken from |Cardelli et al.] ( |1989| ) and 
Gordon et al.| ( |2003| ) respectively. The lines indicate the extinc- 
tion of the diff'erent photometric bands of Table [T] The SMC 
curve only contains values at filter wavelengths explaining the 
strange behavior. Since both J and K bands almost overlap, only 
one of them is indicated in the figure. The extinction in the wave- 
length range used for the SED of J004441.04 is approximately 
the same for both the MW and the SMC. The largest diff'erence 
between both reddening laws are shown in the J and K bands. 
The extinction is small for these IR bands, making that the dif- 
ferent reddening laws hardly diff'er in the wavelength regime of 
interest. Because of the strong resemblance between the MW 
and SMC curve in the used wavelength range, we use the red- 
dening law of the MW to determine the reddening in the SMC. 

The third extinction cause is reddening by the circumstellar 
dust envelope of the post-AGB object itself. We applied a dered- 
dening assuming that the wavelength dependency of the circum- 
stellar extinction is similar to the ISM extinction law of the MW. 
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Fig. 6. Spectrum synthesis of Eu II (top) and Gd II (bottom) lines 
at 6437.640 A and 5733.852 A respectively. The black spectrum 
is J004441.04, the colored spectra are synthetic spectra with dif- 
ferent abundances of the studied element. For more information, 
see text. 



Table 4. SED results of J004441 .04 using a distance of 61 kpc. 



E(B-V) 0.64 + 0.02 
L (Lo) 7600 + 200 
Ur/L^ 0.34 + 0.01 
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Fig. 7. [X/Fe] results of J004441.04, the errorbars represent the 
total uncertainty cr^o?- Some elements are labelled for clarity. For 
more information, see text. 
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Fig. 8. Com parison of [X/Fe] r esults for J004441.04 and 
IRAS06530 ( [Reyniers et aLl|2ob4| ) scaled to [Zr/Fe]. Some el- 
ements are labelled for clarity. The red errorbars of J004441.04 
ratios are the total uncertainty while for IRAS06530, the black 
errorbars represent line-to-line scatter. For more information, see 
text. 



We determined the total reddening by applying a x^ mini- 
malisation on the fit between the scaled model atmosphere and 
the dereddened broadband fluxes. A Monte Carlo simulation is 
used to determine the error on E(B-V) for 1000 arrays with a 
normal distribution of the original flux. For J004441.04, this re- 
sults in an E(B-V) of 0.64 ± 0.02. 



The SED of J004441.04 is shown in Fig. [T0| together with the 
IR spectrum of J004441.04 which was used in |Volk et al.| ( [2QTT ) 
to identify the object as a 21 yum source. Squares indicate the 
dereddened fluxes, triangles indicate original (reddened) fluxes 
and the black line represents the scaled atmospheric Kurucz 
model (see|2]). 

J004441.04 has a clear double-peaked SED where the left 
peak represents the stellar photosphere and the right peak is the 



IR excess created by radiation of dust grains in the ejected AGB 
mass envelope (Fig.[TO]). By integrating the surface of the scaled 
Kurucz model, the luminosity of J004441.04 can be obtained as- 
suming the distance to the SMC is approximately the dist ance to 
J004441.04. We use a distance of 61 kpc ( [Hilditch et aLl|2QQ5] ) 
to determine the luminosity which results in a photospheric lu- 
minosity of 7600 ± 200 L© where the luminosity error is again 
determined via Monte Carlo simulations. The luminosity ratio 
of the IR excess and the photosphere can give some information 
about the strength of the circumstellar reddening. The IR lumi- 
nosity is simply determined by integrating the surface of the IR 
excess. The resulting luminosity ratio Ljr/L^ is 0.34 ± 0.01 so 
about 34% of the radiation emitted by the star is absorbed and 
re-emitted by dust in the circumstellar envelope contributing to 
reddening. We can assume the visible galactic extinction Ay can 
be approximated to Ay ^ 3.1 • E(B-V) and via Fig. |9]we can 



K. De Smedt et al.: Post-AGB stars in the SMC: J00444 1.04-732 136.4 



4.5 
4.0 
3.5 



> 

I 2.5 



:2.0 
1.5 
1.0 
0.5 
0.0 















MW 


_ 




A 


SMC . 


- 






- 






/ 




_ Spitzer 


K 


\ 


1 


1 1 1 



0.0 



0.5 



1.0 1.5 

1/A (l/^m) 



2.0 



2.5 



Fig. 9. Extinction curves of the MW (green) and the SMC (blue). 
The position of the photometric bands used for the SED of 
JQ04 441.04 are indicated. The MW extinction curve is taken 
from [Cardel li et al. ( |1989| ), the SMC extinction curve from 
Gordon et al. (2003). The SMC curve only contains values at 
filter wavelengths explaining the strange behavior. 
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Fig. 11. OGLE lightcurve of J004441 .04 in the I band. For more 
information, see text. 
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Fig. 12. Phase diagram with the lightcurve of J004441.04 (blue 
crosses) folded over the dominant frequency of 97.6 ± 0.3 days 
(red line). 



Fig. 10. SED of J004441.04. Squares indicate the derived dered- 
dened flux, triangles indicate original flux. The black line repre- 
sents the Kurucz atmosphere model used to fit the photometry of 
J004441.04. The red spectrum is the IR spectrum which is used 
in |Volketal.| ( |20TT] ). 



assume the V band extinction of the MW and SMC are approx- 
imately the same. The deduced E(B-V) of 0.64 ± 0.02 makes 
that roughly 20% of the original photosphere flux reaches the 
observer and the other 80% of the photosphere flux is absorbed 
or scattered outside the line-of- sight. Comparing this with the 
luminosity ratio L//?/L^ = 0.34 ± 0.01 indicates that there is a 
significant contribution of the ISM extinction in the SMC and/or 
that there is deviation from spherical symmetry in the circum- 
stellar shell. 

Post-AGB stars pass through the Cepheid II instability strip 
and J004441.04 is no exception. The variability in the I filter 
is displayed in Fig. [TT] A frequency analysis of this light curve 
provides a dominating period of 97.6 ± 0.3 days which is plotted 



in a phase diagram in FigJT2] Blue crosses are the individual 
light curve points of Fig.|ll|and it is clear from Fig. 12 that this 
is indeed the dominant period. Amplitude modulation is clearly 
apparent in Fig.[TT] As the amplitude of the variability is rather 
small (maximum ±0.15 mag in I), the use of stable atmosphere 
models for the spectral analyses is still appropriate. 

We attempted to use the observed period of 97.6 days to con- 
strain the parameters of J004441.04 under the assumpti on that 
this is a pulsation. The linear pulsation cod e of , Wood (19761 
and its updates as described in |Kamath et al.| ( [201Q| ) was used to 
examine the pulsation period of the fundamental mode in stel- 
lar models with luminosities of 7600, 10000, 15000 and 18000 
Lsun and Tgff values between 4700 and 6300 K. The stellar mass 
was obtained from a linear fit between L and M for the two post- 
AGB tracks of Vassiliadis and Wood (1994) for an abundance 
of Z=0.001. The composition of the models was assumed to be 
X=0.75andZ=0.001. 

For our derived Teff of 6250 K for J004441.04, the pulsa- 
tion period of all these models is too small, typically near 50 
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days. This is about half the observed period. A much cooler ef- 
fective temperature is required in order to obtain a model with 
the observed period of 97.6 days. At the luminosity L = 7600 L© 
(our estimated value), a fundamental mode period of 97.6 days 
occurs at Teff ~ 4830 K while at the highest luminosity inves- 
tigated (18000 Lsun) a fundamental mode period of 97.6 days 
occurs at Teff ~ 5010 K. In addition, models with Teff near 6250 
K are pulsationally stable in modes with periods greater than 20 
days while those with Teff near 4900 K are unstable. The pul- 
sation periods are relatively insensitive to composition (X=0.10 
was tried) and the adopted mass. 

The Teff value of ~ 4900K required to get the correct pul- 
sation period is, however, incompatible with the spectrum of 
J004441.04. There are two possible explanations for this dis- 
crepancy. Firstly, it is possible that the pulsation models are 
wrong. These stars are truely extreme: in the inner one third of 
the star radiation provides 90% of the pressure supporting the 
exterior envelope. The stars are difficult to analyse but all the 
relevant physics is included so there should be no obvious prob- 
lem. The metal abundances in the envelope are non-standard so 
that the solar-scaled abundances used will lead to errors in the 
opacity. This is more likely to affect the stability of the models 
than the periods. It is not obvious that the pulsation models have 
significant incomplete physics. 

If the pulsation models give the correct periods, then a sec- 
ond explanation for the discrepancy is required. Here we raise 
the possibility that the observed periodicity is not a pulsation but 
it could be due to other surface phenomena like rapid rotation 
of a spotted star or an orbiting companion. With our adopted Teff 
and L values, and with a stellar mass of 0.65 M©, the breakup pe- 
riod is ~ 94 days, very close to the observed period. So if rotation 
is an explanation for the variability, the star must be rotating at 
near-breakup velocity, which is again incompatible with the de- 
tected broadening of the spectral lines, unless we happen to look 
along the rotation axes. The angle between the line-of- sight and 
the rotation axes should, in that case, be lower than 10°. Also, 
models indicate that a strong rotation and rotational shear in the 
inter shell, may mix ^^N into the ^^C pocket hence impact on the 
efficiency of the s-process nucleosynthesis ( Herwig et al.|2003] 
[Siess et aL|2 004). Considering the other possibility, an orbiting 
companion or dust cloud would have to be very close to the sur- 
face of J004441.04. The SED shows no indication of a hot dust 
excess (Fig. [TO]). All these possibilities do seem very unlikely. 
We are therefore unable to provide a satisfactory explanation for 
the observed periodicity in J004441.04. 



5. Initial mass determination 

Before comparing the results of J004441.04 with theoretical nu- 
cleosynthetic models, its initial mass (mass at the zero-age main- 
sequence) needs to be estimated. Both metallicity and initial 
mass are fundamental parameters for the calculation of models. 
We use t he evolutionary post-AGB tracks of [Vassiliadis & 
|Wood| ( p^94| ) to obtain an accurate estimate of the initial mass. 
The choice of these tracks is based upon the metallicity: these 
authors calculated tracks for a metallicity of Z = 0.001 corre- 
spo nding to the found metallici ty of J004441.04. Unfortunately, 
the [Vassiliadis & Wood] ( |1994| ) tracks start at an eff'ective tem- 
perature of Teff = 10"^ K and hence do not include the 6250 
K of J004441.04. However, the luminosity of post-AGB stars 
remains approximately constant during their transit through the 
HR diagram, allowing linear extrapolation between log Teff and 
log L/L© towards lower temperatures. The need for extrapola- 
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Fig. 13. J004441.04 positioned in the HR diagram based upon 
the results in this paper. The red and green lines respectively rep- 
resent the 1.5 M© and 1.0 M© evolutionary tracks of Vassiliadis 
|& Wood| ( [T994l ) for Z = 0.001. 

tion to lower temperatures may indicate that J004441.04 is in 
the beginning stage of the post-AGB phase. 

Fig. p3] shows the position of J004441.04 in the HR diagram 
together with the 1.5 M© (red line) and 1.0 M © (green line) 
evolutionary tracks of Vassiliadis & Wood ( |1994| ) for Z = 0.001. 
Based upon this figure, we estimate an initial mass of M ^ 1.3 
M© for J004441.04. The tracks for Z=0.001 predict a mass for 
the current post-AGB star of 0.65 M©. With the derived Tg/f and 
L, the radius is 75 R©, and log g = 0.5 also gives M=0.65 M© 
which is consistent with the tracks. 



6. AGB chemical models 

We compare the observed abundance results with predictions 
from two independent stellar evolution codes, with which we 
calculated models of a 1.3M© star of [Fe/H] = -1.4. 



6.1. Mount-Stromlo Evolutionary predictions 

One of the AGB evolutionary models was calculated using the 
same version of the stellar evolution code described in Karakas] 
J2010', and references therein), which uses the [Vassiliadis &| 
[Wood ( 1993j) mass-loss rate on the AGB, and includes the ad- 
dition of C and N-rich lo w-temperature opacities tables from 
Lederer & Aringer ( 2009| ). Convective overshoot is used to in- 
duce the third dredge-up (TDU) in the 1.3M©, Z = 0.0006 
model. The TDU is the inward movement of the convective en- 
velope into regions pr ocessed by pa rtial He-burning during a 
thermal instability (see Herwig 2005 , for a recent review). Some 
form of overshoot is required because low-mass models com- 
puted previously with the same code shows little or no TDU for 
masses < 2 M© at solar metallicity. We calculate one stellar evo- 
lution sequence of a 1.3M© model star with a scaled-solar initial 
abundance pattern with Z = 0.0006. 

We include overshoot by extending the position of the base 
of the convective envelope by 1.0 pressure- scale height. We do 
not vary this free parameter for this study, noting that the final 
predicted carbon abundance of log 6(C) = 8.90 is already 40% 
higher than the observed carbon abundance of log 6(C) = 8.76. 
The final predicted luminosity of the stellar model at the tip of 
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■W- 1.3M3,PMZ=2e-3M,, 
<>- 1.3M^,PMZ=6e-3M,, 
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Fig. 14. The predicted abundance of elements lighter than Fe, 
in log6(X), as a function of atomic number, Z, for a 1.3 M© 
model of Z = 0.0006 ([Fe/H] = -1.4). Abundances are sampled 
at the tip of the AGB, after the last computed thermal pulse. 
Included are the approximate locations (in proton number, Z) 
of some key elements. Predictions are shown for two values of 
the parameter Mmix = 2 x 10"^ (points connected by the solid 
green line) and 6 x 10"^ M© (points connected by the red dashed 
line). This parameter determines the size of the ^^C pocket and ^- 
process enrichment, see text for details. The derived abundance 
of the post- AGB star J004441.04 is shown with error bars. 



7 
6 
5 
4 

X 

<u 3 
o 
2 

1 



-1 



1 

Fe 












1 




r 


1 


— 1 




-¥- 1.3M„PMZ=2e- 


3M, 


». 


















-O- 1.3M,,PMZ=6e- 


3M3 . 


Ni 


















• J004441.04 








- 


|U 


\ll 


Zr 

\ 


!i 




<> 

11 




< 

■ 

i 


fl Dy^^YbHfW 


n < 

r 


1 






1 


r 


M 


1 


f¥|r 


1 1 


1 



30 



40 



50 60 

Z 



70 



80 



Fig. 15. Same as Fig. 14 except for elements heavier than Fe. 



the AGB is 8400 L©, which is achieved after 14 thermal pulses. 
The observed luminosity of 7600 L© is reached after 12 thermal 
pulses. It would be possible to match the observed luminosity 
at the tip of the AGB by either increasing the mass-loss rate or 
increasing the amount of overshoot. 

The ^-process abundance predictions were calculated using 
the post-processing nucle osynthesis code a nd full network of 
320 species described in Lugaro et al."("201 1), with reaction rates 
taken from the JINA REACLIB library ( Cyburt et al._2010 ). For 
the initial c omposition we used t he solar distribution of abun- 
dances from Asplund et al. ( |2009 ) scaled down to [Fe/H] =-1.4. 



Solar abundances of C, N, O, Ne, Mg, S i, S, Ar, and Fe are th e 
pre-solar nebula values from Table 5 of I Asplund et al. ( 2009| ); 
F is the meteoritic value of log 6(F)© = 4.42 from Table 1 
of the same paper (chosen because it has a lower uncertainty), 
and for many of the elements heavier than Fe we use the me- 
teoritic values for the solar abundances (e.g., Sr, Eu, Pb). No 
Qf-enhancement was used at this stage. 

In order to obtain an enrichment of ^-process elements, 
we artificially introduce some protons into the top of the He- 
intershell. Note that this has become standard practice for the 
simple reason that there is not enough ^^C in the H-buming 
ashes of AGB stars to make it an efficient neutron source. Recent 
studies investigating the formation of the ^^C pocket by various 
mechanisms have found that the pro ton abundance in the i nter- 
shell decreases monotonically (e.g., "Cristallo et al.||20()9]). We 



apply the assumption that the proton abundance in the intershell 
decreases monotonically from the envelope value of ^ 0.7 to a 
minimum value of 10""^ at a given point in mass located below 
the base of the envelope. We do this in t he same manner as de- 
scribed in |Alves-Brito et al.| ( |2011| ) and |Kamath et al.| ( |2012| ). 

This partially mixed zone (or PMZ) is required in order to 
facilitate the formation of a ^^C pocket which allows neutrons 
to be released by the^^C(a,n)0^^ reaction ( Straniero et al. 




iGallino et al.|[T998t [Goriely & Mowlavi.,2000, .Herwig. . 

The mass of the proton profile is a free parameter which we set 
as a constant mass. We adopt two choices for the extent of the 
partially mixed proton zone: 1) PMZ = 0.002M© (green solid 
line in Figs.JT4] andJB]), and 2) PMZ = 0.006M© (red dashed 
line in Figs.|14|and|15|). 



In Figs. 14 



and 15 we show the abundances derived for 



J004441.04 in this work along with nucleosynthesis predictions 
from this model. The main point of this comparison is to show 
that a standard, initial scaled-solar abundance pattern used in 
the stellar model produce ^-process overabundances that are not 
compatible with the derived abundances for J004441.04. The 
main shortcomings are that: 

- While the predicted C overabundance is only 40% higher 
than the observed value, the predicted O overabundance is 
too low. This makes that the C/0 prediction is ~18, which is 
clearly in contrast to the detected value of 1.9+0.7. An initial 
alpha enhanced enrichment of [0/Fe] = +0.4, still makes the 
final O abundance too small and the high [0/Fe] = 1.14 indi- 
cates that the star has dredged-up a considerable amount of 
oxygen, above the level found in the stellar model of [0/Fe] 
= 0.35 dex. 

- The total overabundances of the s-process nuclei are clearly 
too low; 

- While the predicted s-process distribution is very similar to 
the observed one, the predicted Pb abundance is significantly 
higher than the observed upper limit (see Sect. [63] and Fig. 

[i3- 

Note also that the final predicted ^^C/^^C ratio is extremely 
high at ^ 1800. An observational constraint on the isotopic ratios 
is not possible with the wavelength coverage we have at present, 
but this should be certainly a priority in our next observational 
season. 

It is plausible that at the low metallicity of [Fe/H] = -1.34, 
J004441.04 evolved from an initially CK-enhanced composition. 
This would explain the discrepancy between the predicted and 
observed Si, S and Ca abundances. 
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6.2. STAREVOL 

To test the uncertainties associated with stellar modeling, an 
additional l.SM©, [Fe/H]=-1.4 model was computed with the 
STAREVOL code (e.g. Siess|20Q7 and references therein). The 
in itial composit i on is slightly different, scaled s olar according 
Cunha ctH] pQQ6!) which is based upon the |Asplund et al. 
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to 



( |2Q05 ) composition with neon enhanced by nearly a factor of 2. g 



This change in the neon abundance was required to reproduce a 
correct seismic solar structure. As before, we do not consider a- 
enhancement and, with our reference composition, the star has 
a metallicity Z=0.0044. We use the Reimers mass loss rate with 
T] = 0.5 up to the beginning of the AGB phase and then switch 
to the |Vassiliadis & Wood| ( |1993] ) prescription. The effect of CO 
opacities is modeled using the analytical fits described in ,Marigo 
( |2Q02| ). Overshooting is also included at the base of the convec- 
tive envelope but no overshoot from the flash-driven convection 
zone into the CO cor e is taken into a ccount. We use the diffu- 
sive approach of Her wig et al.| ( |1999] ) with fover = 0-02. Similar 
to the Stromlo model, this AGB model experiences 15 thermal 
pulses and the TDU starts at pulse number 6. Considering the 
small overshoot parameter, the amount of dredged-up material is 
rather limited and, except for pulse number 9, /I < 0.06 (where A 
is the ratio of the dredged-up mass to the intershell mass increase 
between 2 consecutive pulses). At the tip of the AGB, the stellar 
luminosity is 9400 L©. As previously mentioned, this value can 
be lowered if a higher mass loss rate is used. 

The STAREVOL network includes 55 species and allows 
an accurat e treatment of the nucle osynthesis up to ^'^Cl (for 
details, see |Siess & Amould| ( |2008| ). In this code, mixing and 
nuclear burning are solved simultaneously, once the structure 
has converged. For the s-process nucl eosynthesis, we use the 
post-processing code as described in ,Goriely & Siess| ( [200 1| ) 
which includes about 547 nuclei up to Po with all relevant 
nuclear (n-, p-, a-captures), weak (electron captures, - decays) 
and electromagnetic (photodisintegration) interactions. Nuclear 
reaction rates are taken from the updated Nuclear Astrophysics 
Library of the Brussels University (BRUSLIB, available 
at http://www.astro.ulb.ac.be/bruslib). Following 
[Goriely & M owlavi" (2000), a PMZ of about half the mass 
extension of the pulse (i.e some 6 - 10 x lO'^M©) is introduced 
below the convective envelope at the time of the TDU. In 
this PMZ, the proton profile decreases exponentially from the 
envelope mass fraction down to 10"^ allowing for the formation 
of a substantial ^^C pocket. 

The results of the post-processing calculations are shown by 
the blue line in Fig [16] The model has the same basic shortcom- 
ing, namely an oxygen enrichment ([0/Fe] =0.13) too weak by 
a factor of -10 but a carbon enrichment ([C/Fe] = 1.42) compati- 
ble with the observed value. As illustrated in Fig 16 the surface 
s-process enrichment is also too low with respect to the obser- 
vations. The s-process abundance pattern itself is in good agree- 
ment except again for Pb (see Sect. [63). 

At the end of the AGB evolutionary phase, we are left with a 
relatively small convective envelope of less than 0.05 M©. If we 
assume that our star model experiences a very last thermal pulse 
with a deep TDU, a substantial surface pollution follows and 
a good agreement with the observations can be achieved (red- 
solid line in Fig. 18). This last-pulse modeling corresponds to a 
post-processing calculation identical to those characterizing the 
previous TDU episodes but for which the TDU is just assumed 
to extend to about 7 x 10"^ M© deep into the C-rich region, 
i.e about half of the pulse mass is now diluted into the enve- 
lope. This improved fit to the observation is due to the small 
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-■- AGB model 
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Fig. 16. Predicted elemental distribution in log e (top) and [X/Fe] 
(bottom panel) as a function of the atomic number Z for a 
1.3M0,[Fe/H] = -1.4 model computed with STAREVOL. In each 
panel, the blue-dashed line corresponds to the post-processed 
AGB model including 15 thermal pulses and the red line shows 
the results expected from the occurrence of an additional final 
thermal pulse followed by a deep TDU. 



remaining envelope mass (M^nv < 0.05Mo) which limits the 
dilution of the dredged-up material. The observed distribution 
of s-process elements can be nicely reproduced except for Pb 
which, in both simulations, is overproduced (see Sect. 6.3 ), as 
classically done in all s-pr ocess simulations at low-metallicity 
( [Goriely & Mowlavi|2000| ). Concerning the light elements, this 
last dredge-up brings significantly more C, O, F, Ne and Na to 
the surface. The C is now in disagreement with the observations 
but on the other hand, we can reproduce the oxygen enrichment, 
yielding a similar C/0~ 20 as before. We also note that the un- 
derabundance of the a-elements Si, S and Mn remains. It is inter- 
esting to see how the use of diff'erent abundances units, namely 
log(6) and [X/Fe] can be misleading. Finally, note that within the 
s-process model adopted here, the partial mixing of protons into 
the C-rich region also leads to a significant production of F, Ne 
and Na (see e.g [Goriely & Mowlavi| ( [2000 )) which is also the 
case for the Mount- Stromlo models as shown in Fig. 18 The es- 
timated abundance of Na by the models is seen in Fig. 18 to be in 
rather good agreement with observation. It would be extremely 
valuable to confirm observationally the predicted F overproduc- 
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Fig. 17. Pb I line region in the spectrum of J004441.04 com- 
pared with model predictions. In green, the Mount Stromlo 
model prediction (Section 6.1) and in red and blue the predic- 
tions obtained by the STAREVOL code (Section 6.2). 

tion of about 2 dex. Unfortunately, at these photospheric tem- 
peratures, F has no suitable lines in the sampled spectrum and 
its abundance can not be quantified. 

6.3. Pb discrepancy 

Figs. [15] [16] and [18] clearly show that current stellar evolution 
codes predict strong Pb overabundances for J004441.04. The 
strongest predicted Pb line in our sampled spectral range is 
4057.807 A. Unfortunately, this spectral region has a low S/N 
preventing an accurate Pb abundance determination. We focused 
on this spectral region and looked carefully on the three expo- 
sures. We obtained a weighted average of the three UVES spec- 
tra and then confronted this final spectrum with spectral synthe- 
sis models. We used the Pb abundances predicted by the chem- 
ical models in Fig. [TT] The black spectrum is the normalised 
weighted mean spectrum of J004441.04, the green spectrum dis- 
plays the predicted Pb abundance of the Mount- Stromlo model 
with PMZ = O.OO6M0, the red and blue spectrum represent the 
STAREVOL predictions with the blue spectrum the model with 
a last thermal pulse. Although the stellar spectrum is of poor 
quality, it is clear from Fig. 17 that the predicted strong line is 
not present in the stellar spectrum. From Fig. [TT] we estimate 
an upper Pb abundance of 3.00 for J004441.04. Unfortunately 
we do not have a spectrum of IRAS06530 in this spectral do- 
main. In the late thermal pulse scenario (Fig. j_8), the predicted 
Pb abundance becomes very large, which is incompatible with 
the detected spectrum. To reduce the production of Pb, the neu- 
tron in the PMZ production must be reduced. Partial pollution of 
the ^^C pocket by ^"^N induced by rotational mixing (e.g. Herwig 
|eFa l. 2003 ; Siess et al.|2004| ) may be needed to solve this prob- 
lem but this will need a detailed investigation which is outside 
the scope of this paper. 

7. Discussion and Conclusion 

We presented a detailed analysis of J004441.04, the only known 
21 yum object of the SMC. This metal poor star ([Fe/H] = -1.34 
± 0.32), turns out to be among the most s-process enriched stars 
known to date while displaying only a moderate C/0 ratio of 1.9 



± 0.7. The atmospheric parameters (Te/f = 6250 ± 250, log g 
= 0.5 ± 0.5) combined with the luminosity of 7600 ± 200 L© 
as well as the pulsation period of 97.6 ± 0.3 days, show that we 
witness the post-AGB phase of a star of low initial mass. This 
initial mass is estimated to be ^ 1.3 M©. Our findings on an extra- 
Galactic source corroborate the conclusion that 21 yum stars are 
post-Carbon stars. Although of lower metallicity, the object dis- 
plays very similar overabundances as the most s-process rich 
post-AGB star known in our Galaxy (IRAS06530-0213). 

Our AGE model predictions are based on model calculations 
which include a forced overshoot of the convective boundary as 
well as an artificial inclusion of a proton profile into the inter- 
shell. We used two independent codes for these predictions. For 
comparsion, we also included in Fig. 18 the nearest model in the 
coarser grid of models by [Cristallo et al.| (201 1 ) (1.5 M©, Z = 
0.001) which yields similar discrepancy like our two models. 

While the predicted C overabundance is compatible with the 
observations, the predicted O abundance is significantly lower 
resulting in a predicted C/0 ratio of ~20 which is clearly too 
high. More oxygen could be dredged-up in the models, if one 
includes additional over shooting below the thermal pulse, as 
shown by |Herwig| ( [2Q00] ). 

While the models fit the s-process distribution well, as de- 
picted in Fig. 
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the absolute abundances are not well matched. 
In our late thermal pulse scenario, we assume that a thermal 
pulse occurred with dredge-up and a limited dillution by the re- 
maining reduced envelope mass. In this model, the very high 
overabundances are matched better. A noticeable exception is 
that the predicted Pb abundance is much larger than the detected 
upper limit (Fig.[T7]). A higher S/N spectrum in the region of Pb 
lines is needed to better quantify the Pb abundance which would 
yield an ever stronger test on the nucleosynthesis but our spec- 
trum is incompatible with such high Pb abundance. With the data 
at hand, no isotopic ratio could be determined yet but the pre- 
dicted ^^C/^^C of ~ 1800 is extremely high. The nearest model 
in the coarser grid of models by |Cristallo et al.| ( [201 Ij yields 
similar discrepancy with even a higher predicted C/0 ratio and 
i^C/i^C values. 

Our calculations demonstrate that there is only a weak de- 
pendency of the theoretical predictions on the adopted stellar 
evolution code ! 

It is interesting to compare the abundance pattern of 
J004441.04 with the patterns observed in peculiar metal poor 
stars like the CEMP-s/r stars, that is, a carbon-enhanced metal- 
poor star with en richments in both s and r-process elements 
(see discussi ons in [Beers & Christlieb|2005"| [Jonsell et al.|2006| 



[Sneden et ar]|2008[ [Lugaro et al.||2011| ). According to the def- 
inition given in Jonsell et al. (2006), CEMP-s/r stars have 
[Eu/Fe]>l and [Ba/Eu]>0. The post-AGB star J004441.04 sat- 
isfies all of these criteria. It is carbon rich with [C/Fe] = 1.67 
and [Ba/Eu] > 0, using the [Ba/Fe] = 2.61 estimated earlier. The 
estimated Ba (or observed La, which is a good proxy for ele- 
ments at the second ^-process peak) and Eu abundances put this 
star well into the region of the [Eu/Fe] -[ Ba/Fe] plane occupied 
by the CEMP-s/r stars (see Fig. 5 from Lugaro et al. 201 1| ). 
Abundance wise, J004441.04 could be seen as a more metal rich 
analogue of the CEMP-s/r stars. However, unlike many of the 
stars within the CEMP-s/r class, the distribution of heavy ele- 
ments in J004441.04 is similar to that of a pure s-process, except 
for the anomalously low Pb abundance. What is unusual about 
J004441.04 is how the star obtained such large enrichments of 
heavy elements, while keeping the C/0 ratio low. 

It is fair to say that many of the chemical, physical and 
pulsational properties of the interesting SMC post-AGB star 
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Fig. 18. The abundance pattern of J004441.04 in compari- 
son with the model predictions. In red, the Mount Stromlo 
model (Section 6.1), in blue the predictions obtained by the 
STAREVOL code (Section 6.2), the green dashed curve is com- 
puting assuming a late thermal pulse with deep dredge-up occurs 
when the diluti on is small. The mage nta line is the 1.5 M©, Z = 
0.001 model of |Cristallo et al.| ( |2QTT| . 



J004441.04 are, as yet, unaccounted for by our model predic- 
tions. This study is the first of a series in which we aim at a 
systematic study of post- AGB stars in the Magellanic clouds. 
These objects with constrained distances will provide unprece- 
dented and systematic tests of the AGB model predictions, com- 
plementary to the limited stu dy of AGB stars (e.g. |de Laveriiy 
|et al.||2006| |Abia et aL]|2011| ) for which much lower overabun- 
dances are detected. The ultimate aim of this programme is to 
make progress in the understanding of the complex interplay be- 
tween mixing, nucleosynthesis and mass-loss which characterise 
the final evolution of solar-like stars. 



Acknowledgements. HVW and KDS acknowledge financial support from the 
Research Council of K.U.Leuven under grant GOA/2008/04 and from the 
Scientific Fund of Flanders (FWO) under the grants G.0703.08 and G.0470.07. 
AIK is grateful for the support of the NCI National FaciHty at the ANU. AIK 
is an ARC Future Fellow and is supported under grant FTl 01 00475. LS and SG 
are FNRS research associates. The authors thank Chris Sneden for his comments 
on the manuscript. 



References 

Abia, C, Cunha, K., Cristallo, S., et al. 2011, ApJ, 737, L8 

Abia, C, de Laverny, R, & Wahlin, R. 2008, A&A, 481, 161 

Abia, C, Dominguez, L, Gallino, R., et al. 2002, ApJ, 579, 817 

Alves-Brito, A., Karakas, A. I., Yong, D., Melendez, J., & Vasquez, S. 2011, 

A&A, 536, A40 
Angelou, G. C, Church, R. R, StancHff'e, R. J., Lattanzio, J. C, & Smith, G. H. 

2011, ApJ, 728, 79 
Asplund, M., Grevesse, N., & Sauval, A. J. 2005, in Astronomical Society of 

the Pacific Conference Series, Vol. 336, Cosmic Abundances as Records of 

Stellar Evolution and Nucleosynthesis, ed. T. G. Barnes III & F. N. Bash, 25 
Asplund, M., Grevesse, N., Sauval, A. J., & Scott, R 2009, ARA&A, 47, 481 
Beers, T. C. & Christlieb, N. 2005, ARA&A, 43, 531 
Bolatto, A. D., Simon, J. D., Stanimirovic, S., et al. 2007, ApJ, 655, 212 
Busso, M., Palmerini, S., Maiorca, E., et al. 2010, ApJ, 717, L47 
Cardelli, J. A., Clayton, G. C, & Mathis, J. S. 1989, ApJ, 345, 245 
Castelli, F & Kurucz, R. L. 2004, ArXiv Astro-ph/0405087 
Church, R. R, Cristallo, S., Lattanzio, J. C, et al. 2009, PASA, 26, 217 
Coelho, R, Barbuy, B., Melendez, J., Schiavon, R. R, & Castilho, B. V. 2005, 

A&A, 443, 735 
Cristallo, S., Piersanti, L., Straniero, O., et al. 2011, ArXiv e-prints 



Cristallo, S., Straniero, O., Gallino, R., et al. 2009, ApJ, 696, 797 

Cunha, K., Hubeny, I., & Lanz, T. 2006, ApJ, 647, L143 

Cutri, R. M., Skrutskie, M. R, van Dyk, S., et al. 2003, VizieR Online Data 

Catalog, 2246, 
Cyburt, R. H., Amthor, A. M., Ferguson, R., et al. 2010, ApJS, 189, 240 
de Laverny, R, Abia, C, Dominguez, I., et al. 2006, A&A, 446, 1107 
Dekker, H., D'Odorico, S., Kaufer, A., Delabre, B., & Kotzlowski, H. 2000, 

in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference 

Series, Vol. 4008, Society of Photo-Optical Instrumentation Engineers (SPIE) 

Conference Series, ed. M. lye & A. F. Moorwood, 534-545 
Fouque, R, Chevallier, L., Cohen, M., et al. 2000, A&AS, 141, 313 
Gallino, R., Arlandini, C, Busso, M., et al. 1998, ApJ, 497, 388 
Giridhar, S., Ferro, A., & Parrao, L. 1997, PASP, 109, 1077 
Gordon, K. D., Clayton, G. C, Misselt, K. A., Landolt, A. U., & Wolff', M. J. 

2003, ApJ, 594, 279 
Goriely, S. & Mowlavi, N. 2000, A&A, 362, 599 
Goriely, S. & Siess, L. 2001, A&A, 378, L25 

Habing, H. J. & Olofsson, H., eds. 2003, Asymptotic giant branch stars 
Herwig, F 2000, A&A, 360, 952 
Herwig, F 2005, ARA&A, 43, 435 

Herwig, F, Blocker, T., Langer, N., & Driebe, T. 1999, A&A, 349, L5 
Herwig, F, Langer, N., & Lugaro, M. 2003, ApJ, 593, 1056 
Hilditch, R. W., Howarth, I. D., & Harries, T. J. 2005, MNRAS, 357, 304 
Hrivnak, B. J., Volk, K., Geballe, T. R., & Kwok, S. 2008, in lAU Symposium, 

Vol. 251, lAU Symposium, ed. S. Kwok & S. Sandford, 213-214 
Hrivnak, B. J., Volk, K., & Kwok, S. 2009, ApJ, 694, 1147 
Jonsell, K., Barklem, R S., Gustafsson, B., et al. 2006, A&A, 451, 651 
Kamath, D., Karakas, A. I., & Wood, R R. 2012, ApJ, 746, 20 
Kamath, D., Wood, R, Van Winckel, H., & van Aarle, E. 201 1, in 9th Pacific Rim 

Conference on Stellar Astrophysics, ed. , ASP conference series, in press 
Kamath, D., Wood, R R., Soszyhski, I., & Lebzelter, T. 2010, MNRAS, 408, 522 
Karakas, A. I. 2010, MNRAS, 403, 1413 

Karakas, A. I., Campbell, S. W., & Stancliff'e, R. J. 2010, ApJ, 713, 374 
Kiselman, D. 2002, Highlights of Astronomy, 12, 429 
Kobayashi, C, Karakas, A. I., & Umeda, H. 2011, MNRAS, 414, 3231 
Kupka, F., Piskunov, N., Ryabchikova, T. A., Stempels, H. C, & Weiss, W. W 

1999,A&AS, 138, 119 
Kwok, S., Volk, K. M., & Hrivnak, B. J. 1989, ApJ, 345, L51 
Lambert, D. L., Luck, R. E., & Bond, H. E. 1983, PASP, 95, 413 
Lederer, M. T. & Aringer, B. 2009, A&A, 494, 403 
Lobel, A. 2006, in lAU Joint Discussion, Vol. 4, lAU Joint Discussion 
Luck, R. E., Moff'ett, T. J., Barnes, III, T. G., & Gieren, W. R 1998, AJ, 115, 605 
Lugaro, M., Karakas, A. I., Stancliff'e, R. J., & Rijs, C. 2011, ArXiv e-prints 
Marigo, R 2002, A&A, 387, 507 

Posch, T., Mutschke, H., & Andersen, A. 2004, ApJ, 616, 1167 
Reddy, B. E., Lambert, D. L., Gonzalez, G., & Yong, D. 2002, ApJ, 564, 482 
Reyniers, M. & Van Winckel, H. 2003, A&A, 408, L33 
Reyniers, M., Van Winckel, H., Gallino, R., & Straniero, O. 2004, A&A, 417, 

269 
Richter, O.-G., Tammann, G. A., & Huchtmeier, W K. 1987, A&A, 171, 33 
Romano, D., Karakas, A. I., Tosi, M., & Matteucci, F. 2010, A&A, 522, A32 
Schlegel, D. J., Finkbeiner, D. R, & Davis, M. 1998, ApJ, 500, 525 
Siess, L. 2007, A&A, 476, 893 
Siess, L. & Arnould, M. 2008, A&A, 489, 395 
Siess, L., Goriely, S., & Langer, N. 2003, PASA, 20, 371 
Siess, L., Goriely, S., & Langer, N. 2004, A&A, 415, 1089 
Sneden, C, Cowan, J. J., & Gallino, R. 2008, ARA&A, 46, 241 
Sneden, C. A. 1973, PhD thesis. The Univerisity of Texas at Austin. 
Stancliff'e, R. J., Glebbeek, E., Izzard, R. G., & Pols, O. R. 2007, A&A, 464, L57 
Straniero, O., Gallino, R., Busso, M., et al. 1995, ApJ, 440, L85 
Szymanski, M. K. 2005, Acta Astron., 55, 43 

Udalski, A., Kubiak, M., & Szymanski, M. 1997, Acta Astron., 47, 319 
van Aarle, E., van Winckel, H., Lloyd Evans, T., et al. 201 1, A&A, 530, A90 
Van Winckel, H. 2003, ARA&A, 41, 391 
Van Winckel, H. & Reyniers, M. 2000, A&A, 354, 135 
Vassiliadis, E. & Wood, R R. 1993, ApJ, 413, 641 
Vassiliadis, E. & Wood, R R. 1994, ApJS, 92, 125 
Volk, K., Hrivnak, B. J., Matsuura, M., et al. 2011, ApJ, 735, 127 
Wood, R R. 1976, MNRAS, 174, 531 
Zaritsky, D., Harris, J., Thompson, I. B., Grebel, E. K., & Massey, R 2002, AJ, 

123, 855 



14 K. De Smedt et al.: Post-AGB stars in the SMC: J00444 1.04-732 136.4 

Table 5. Linelist used for the abundance determination of J004441.04. The right 
column shows the calculated equivalent width. EWs indicated with SS are calcu- 
lated via synthetic modelling. 



ion 


Z 


wavelength (A) 


EP(eV) 


loggf 


EW (mA) 


CI 


6 


6655.517 


8.54 


-1.94 


56.0 






6012.225 


8.64 


-2.00 


48.8 






6014.834 


8.64 


-1.59 


71.0 






6016.444 


8.64 


-1.83 


74.6 






5551.578 


8.64 


-1.90 


53.2 






6002.987 


8.65 


-2.17 


39.3 


01 


8 


6155.971 


10.74 


-0.67 


18.7 






6454.444 


10.74 


-1.08 


9.7 






6453.640 


10.74 


-1.30 


7.3 


Nal 


11 


6154.230 


2.10 


-1.56 


11.4 






5688.210 


2.10 


-0.42 


74.7 


Mgl 


12 


5711.100 


4.34 


-1.75 


38.0 (SS) 


Si I 


14 


6237.330 


5.61 


-0.53 


35.6 






5948.550 


5.08 


-1.22 


22.5 


SI 


16 


6757.160 


7.87 


-0.20 


20.7 






6748.790 


7.87 


-0.35 


10.5 


Cal 


20 


6102.727 


1.88 


-0.80 


59.7 






6449.820 


2.52 


-0.62 


33.3 






5590.130 


2.52 


-0.74 


26.6 






6439.075 


2.53 


0.39 


116.4 


Sell 


21 


6604.600 


1.36 


-1.23 


98.6 






5552.240 


1.45 


-2.08 


15.3 






5667.153 


1.50 


-1.21 


68.3 


Till 


22 


4798.521 


1.08 


-2.67 


84.5 






6606.980 


2.06 


-2.85 


29.0 






5010.212 


3.10 


-1.30 


76.7 


CrII 


24 


5510.730 


3.83 


-2.48 


18.2 






5310.700 


4.07 


-2.28 


15.8 






5502.090 


4.17 


-1.99 


22.3 


Mnl 


25 


6021.800 


3.07 


0.03 


39.6 






6016.650 


3.07 


-0.22 


33.1 


Fel 


26 


5198.711 


2.22 


-2.14 


33.3 






5049.820 


2.28 


-1.35 


80.7 






6421.351 


2.28 


-2.01 


33.4 






6252.555 


2.40 


-1.72 


41.7 






6191.558 


2.43 


-1.42 


60.1 






6136.596 


2.45 


-1.40 


65.5 






6230.723 


2.56 


-1.28 


59.2 






6137.692 


2.59 


-1.40 


47.5 






5281.790 


3.04 


-0.83 


76.9 






5324.179 


3.24 


-0.10 


108.5 






5615.644 


3.33 


0.05 


118.6 






6400.001 


3.60 


-0.29 


66.2 






6411.649 


3.65 


-0.66 


37.7 






5445.042 


4.39 


0.04 


40.0 






6024.058 


4.55 


-0.06 


29.5 






6419.950 


4.73 


-0.09 


15.3 






5572.842 


3.40 


-0.28 


91.9 


Fell 


26 


5991.368 


3.15 


-3.56 


54.7 






6084.099 


3.20 


-3.80 


48.5 






6147.741 


3.89 


-2.73 


66.1 






6383.721 


5.55 


-2.14 


17.6 


Nil 


28 


6767.780 


1.83 


-2.17 


11.6 






5035.370 


3.63 


0.29 


35.1 






6772.320 


3.66 


-0.98 


5.2 






5694.990 


4.09 


-0.61 


5.1 


Znl 


30 


4810.540 


4.08 


-0.17 


72.5 (SS) 


YII 


39 


5289.815 


1.03 


-1.85 


144.7 
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Table 5. continued. 



ion 


Z 


wavelength (A) 


EP (eV) 


loggf 


EW (mA) 






5728.890 


1.84 


-1.12 


149.4 


Zrll 


40 


5124.982 


1.53 


-1.50 


117.2 






5477.822 


1.83 


-1.40 


86.2 


La II 


57 


5482.268 


0.00 


-2.06 


158.9 






5712.391 


0.17 


-1.96 


147.7 






5936.210 


0.17 


-2.07 


118.5 






5880.633 


0.24 


-1.92 


158.3 






6146.523 


0.24 


-2.47 


83.2 






5062.918 


0.77 


-1.72 


91.9 






6126.075 


1.25 


-1.19 


120.6 






5566.925 


2.38 


-1.00 


51.1 






6399.030 


2.65 


-0.52 


70.2 


Cell 


58 


5031.986 


1.41 


-0.90 


62.5 






5613.694 


1.42 


-1.00 


67.6 






5459.193 


1.62 


-0.58 


94.4 






5959.688 


1.63 


-0.84 


47.6 






6143.376 


1.70 


-0.80 


70.0 






6098.326 


1.77 


-0.61 


61.3 






5359.508 


1.78 


-0.78 


60.1 






5685.836 


1.90 


-0.43 


54.4 


PrII 


59 


5135.140 


0.95 


-0.13 


118.6 






5605.642 


0.96 


-0.65 


71.3 


Ndll 


60 


5221.572 


0.38 


-1.34 


120.7 






5934.738 


0.75 


-1.39 


109.7 






6365.540 


0.93 


-1.37 


77.2 






6031.270 


1.28 


-0.83 


103.8 






6803.980 


1.44 


-0.79 


108.2 






6737.763 


1.60 


-0.78 


93.4 






6650.517 


1.95 


-0.32 


105.2 


Smll 


62 


4791.580 


0.10 


-1.24 


111.2 






6731.813 


1.17 


-0.52 


88.0 


EuII 


63 


6437.640 


1.32 


-0.27 


87.2 (SS) 


Gdll 


64 


5733.852 


1.37 


-0.89 


26.3 (SS) 


Dyll 


66 


5169.688 


0.10 


-1.95 


68.3 (SS) 


ErII 


68 


5229.319 


2.60 


-0.39 


16.5 (SS) 


Ybll 


70 


5352.954 


3.75 


-0.34 


53.6 (SS) 


Lull 


71 


6463.107 


1.46 


-1.05 


126.0 (SS) 


Hfll 


72 


4790.708 


2.20 


-1.11 


38.5 (SS) 


WII 


74 


5104.432 


2.36 


-0.91 


46.3 (SS) 



